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We investigate the influence of an external magnetic field on chiral symmetry breaking 
in a four-fermion interaction model at finite temperature and chemical potential. Using the 
Fock-Schwinger proper time method, we calculate the effective potential for the four-fermion 
' interaction model to leading order in the 1/N C expansion. The phase structure of the chiral 

symmetry breaking is determined in the T-^t, H-T and [i-H planes. The external magnetic 
field modifies the phase structure. It is found that a new phase appears for a large chemical 
■ potential. 
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Q . §1. Introduction 

O , 

The strong interaction between quarks and gluons is described by QCD. In recent 
^r) . years, much interest has been paid to the phase structure of QCD. The Lagrangian of 

QCD has a SU (N) l <8> SU(N)r symmetry for massless quarks. The chiral SU (N) l ® 
SU(N)r symmetry is broken down to the flavor SU(N) symmetry by the QCD 
Ph, dynamics. It is expected that the chiral symmetry is restored at high temperature 

and /or high density. A phase transition takes place from hadronic matter to a quark- 
gluon plasma or a color super-conducting phase under hot and dense conditions. This 
affects the structure of spacetime and the evolution of stars. 

A state of high temperature and density was realized in the early states of our 
^ ■ universe. There is a possibility that the early universe contained a large primordial 

magnetic field. Another dense state at low temperature is found in the cores of 
neutron stars. Some neutron stars, called pulsars, possess a strong magnetic field 
and emit high energy radiation. The existence of a so-called 'magnetar', which is 
a star with a much stronger magnetic field, has been reported. To describe such 
state, it is very important to investigate thermal and magnetic effects on the QCD 
ground state. To study dense states, we introduce a chemical potential. The phase 
structure of QCD is studied in an external magnetic field H at finite temperature T 
and chemical potential fx. 

Because a long-range correlation plans an essential role in the phase transition, 
the QCD ground state is determined non-perturbatively. There are many methods 
to deal with non-perturbative aspects of QCD, e.g. lattice QCD, the Schwinger- 
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Dyson equation, the strong coupling expansion, and so on. In the present paper, we 
consider a four-fermion interaction model as a low-energy effective theory of QCD. 1 ) 
It is known that this model has properties similar to those of non-perturbative QCD 
at the low energy regime. In the model, chiral symmetry is dynamically broken down 
through fermion and anti-fermion condensation for a sufficiently strong four-fermion 
interaction. 

Employing the Matsubara formalism of finite temperature field theory, the ther- 
mal effect on chiral symmetry breaking was investigated in the Nambu-Jona-Lasinio 
(NJL) model. 2 ) For a high temperature and/or a large chemical potential, the expec- 
tation value of a composite state consisting of a fermion and an anti-fermion, (•tp^p), 
disappears, and broken chiral symmetry is restored. A tricritical point between the 
first and second order phase transitions was observed in the critical curve in the T-fi 
plane for H = 0. 3 )~ 5 ) 

A method to treat a strong external magnetic field was developed by Schwinger 
in 1951. 6 ) This method is called the Fock-Schwinger proper time method. In Refs. 
7) and 8), the proper time method is applied to a four-fermion interaction model 
to study the influence of a magnetic field. It is shown that an external magnetic 
field has the effect of effectively reducing the number of spacetime dimensions and 
enhancing the chiral symmetry. Thus, the expectation value (ipip) slightly increases 
in a magnetic field at T = /j, = 0. This method has been extended to a variety of 
situations with external fields. 9 )~ 12 ) 

Here we apply the Fock-Schwinger proper time method to the NJL model at 
finite temperature, chemical potential and external magnetic field. Using an explicit 
expression of the fermion Green function in an external magnetic field at finite T 
and //, 13 ) we calculate the effective potential to leading order in the l/iV c -expansion. 
We use the observed quantities and to fix the parameters in the model. The 
ground state is determined by observing the minimum of the effective potential. 
Evaluating it numerically, we determine the phase structure of the theory in the 
T-/x, H-T and [i-H planes and investigate the combined effects of the magnetic field, 
temperature and chemical potential. Finally, we give some concluding remarks. 

§2. Effective potential with T, /x and H 

In this section, we briefly review the fermion Green function at finite T and 
fi with an external magnetic field and calculate the effective potential for a four- 
fermion interaction model to leading order in the 1/N C expansion. There are some 
works 3 )' 14 )' 17 )' 18 ) in which the original Fock-Schwinger proper time method 6 )' 15 ) is 
applied to finite temperature and/or chemical potential situations. It has been 
pointed out that a naive Wick rotation is not always valid for a proper time in- 
tegral. In order to investigate physical situations, we carefully choose the proper 
time contour, as demonstrated in Ref. 13). We apply this result to the calculation 
of the effective potential in the four-fermion interaction model. 

Here we consider the NJL model, which is well known as one of the low energy 
effective theories, in the chiral limit of QCD. The Lagrangian of the NJL model is 
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given by 

^2 i>i^i J + j ^2 tpiijsTijip. 

>i=u,d J \i,j=u,d 



i=u,d 




where N c is the number of colors, the indices i and j denote the fermion flavors, 
Qi represents the electric charge of the quark fields (Q u = 2e/3, Qd = — e/3), and 
G is the effective coupling constant. We consider two flavors of quark fields. The 
quark fields ipi belong to the fundamental representation of the color SU(N C ) group 
and a flavor isodoublet. f represents the isospin Pauli matrices. This Lagrangian 

is invariant under the global chiral transformation ipi — > exp|i6»^7 5 | Be- 
low, we do not include the flavor indices explicitly. In practical calculations, it is 
more convenient to introduce the auxiliary field a and tt, and hence we consider the 
Lagrangian 

C = i> (ifl - QA ~ o ~ H5T ■ tt) V + ^ + vf 2 ) . (2) 

From the equations of motion for a and tt, we obtain the correspondences a ~ 
-(G/N c )^ and vr ~ -(G/N c )ipij 5 fip. 

To find the ground state of the NJL model at finite T, and H, we evaluate 
the effective potential. Following the standard procedure of the imaginary time 
formalism, 16 ^ we introduce the temperature and the chemical potential. Here, we 
assume that the ground state does not break the local £7(1) em symmetry, i.e. (tt + ) = 
(ir~) = 0. Under this assumption, we can choose the isospin singlet ground state 

from the degenerate state using the chiral transformation ip — > exp |i#-7j-75 1 V 7 - To 

leading order in the 1/N C expansion, the effective potential V e g(a,TT = 0) is given by 

1 2 1 m 1 ( if) + Qfi — ^74 — a 



V„( g , f - 0) - -„> - — TVh [ - j , (3) 

where we assume fj, = fj, u = fj.^ and Tr denotes the trace over spacetime coordinates, 
spinor and flavor indices. The factor (3V in the denominator comes from the 4- 
dimensional volume, 

ff3 poo 

dxi / d 3 x = (3V, (4) 



J-oo 



where j3 = 1/(&bT), with ks the Boltzmann constant. Below, we set ks = 1. The 
effective potential © is normalized to satisfy V e g(cr = 0) = 0. 

The second term on the right-hand side of Eq. © is rewritten using the relation 



Trln( \^ + Q4-^7 4 CJ )=Tr/ < SV "'- ; '" i ' / "'- (5 > 



o 



*' Because the electric charges for up and down quarks are different, the flavor SU(N) symmetry 
is broken down explicitly. Only the symmetry under this transformation remains. 
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where S(x, y; m) is the fermion Green function, which satisfies the Dirac equation 

(ifi + Qfi- ^74 - m) S(x,y;m) = 5%{x - y) , (6) 

where 5^(x — y) = 5(x4 — y4)5 3 (x — y). 

To solve the Dirac equation ©, we introduce a bispinor function G(x,y;m), 
defined by 

S(x, y; m) = {if) + - ^74 + m) G(x, y; m) . (7) 
Substituting Eq. into Eq. ©, we obtain 



W^.dj ~ iQAjf + \QF jk oi k - (id 4 - iiif - m 2 J G(x, y; 



;m)=6£(x-y), (8) 



where Fij is the three-dimensional field strength and a^ k = (i/2)[ r y : > , 7*]. 
Performing the Fourier series expansion, 



G{x,y;m) = - ]T e -^^^)G n (f,y» , co n = ( 2n + 1 ) 7r , (9) 

n=— 00 

we rewrite Eq. (|SJ) into the form 
"• .2 { l u . 

172 - 727iJ 



-(u n -in) -m >G n (x,y;m) = S (x-y), (10) 

where we consider a constant magnetic field along the a^-direction, i.e. F12 = —F21 = 
H, for simplicity.*) 

As is well known, the solution of Eq. 1)101) can be obtained using the Fock- 
Schwinger proper time method. 6 )' 15 ) In Ref. 14), it is pointed out that the original 
Fock-Schwinger method must be modified in the case that T / and \i 7^ 0. The 
function G n (x,y;m) is determined by 

~ J e" 5 ^/ 4 fj° U n (x,y;r) dr for n > 

G n (x, y; m) - j ^ i/A ^ ^ ^ ^ ^ ^ < Q , (11) 

where U n {r) is the proper time evolution operator derived to be 

= (4ff )3/2| T |3/2 Xr) ^ {* Q X 
i 

x exp 



x - y)iQFij[coth(QFr)]jk(x - y) k 

-it 1 -QFjkcrjk - K - ^) 2 - m2 



(12) 



*' If we consider a constant electric field, we obtain an equation similar to Eq. 11UL except for a 
independent potential. Due to the presence of this z-dependent term, we cannot analytically solve 
this equation in an external electric field at finite T and fi. 



J_ 
A 2 



Im[f ] 



m7l 



(m=l,2,3 ; ) 




c, 



Re[f ] 



Fig. 1. The contours of the integration in Eq. 113|l . The circles denote the poles that come from 
r~ 3 / 2 and 1/ sin (QHr). A is a proper time cut-off. 



Substituting Eqs. © and (|llj) into Eq. (J7J), we obtain the fermion Green function 
S(x,y;m). It should be noted that there is a problem of ambiguity in selecting the 
proper time contour for n > T > 0. 17 ^ We must be careful to take the physical 
contour. 

Performing the integration over m in Eq. ©, we obtain the explicit form of the 
effective potential 



1 „-37ri/4 °° 



where (c.c.) denotes the complex conjugate of the second term. In the limit \i — > 0, 
the effective potential (j!3fl reproduces the previous result obtained in Ref. 18). In 
the other limit, T — > 0, Eq. ()13j) is consistent with the results given in Refs. 4) and 
19). 

The integrand in Eq. (|13|) has many poles, which come from | -7~ | 3 / 2 and cot (QHt) 
on the real axis. The contour is determined by the boundary conditions of the Green 
function. The physical contour is obtained in Ref. 6) for T = \i = 0. Because our 
result must coincide with that obtained using this contour in the simultaneous limits 
T — > and fj, — * 0, the contour for the proper time integral in the second term of 
Eq. lfT5J) should be C\ in Fig. ^ We should employ the contour C2 for the complex 
conjugate term. 

To determine the contribution from the poles on the real r axis, we rewrite Eq. 
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(|13|) in the form 



a 2 



V&i?) = ^ + 1(a) - I(a = 0) + I* (a) -I*(a = 0), (14) 
where 

1(a) = (^E/ Ci ^^ C ot(Q^)e--(--) 2 — 2 . (15) 

After the Wick rotation in the complex T-plane, the contribution from the poles is 
classified into the following two cases. 

a) uj 2 + a 2 - fi 2 = (irT) 2 + a 2 - n 2 > 

In this case, the integrand in Eq. I)15|) converges as r — ► —too. Hence, the 
expression (jl5|) can be continued to imaginary r without encountering any poles: 

1 °° /"OO 

I(a)+I*(a) = - TJ -Y. drf(a-n;r), (16) 



f(a; n;r) = ^ coth(Q^r) co S (2uo n ^r)e- T ^ 2 ~^ + ^ . (17) 

We can rewrite the summation in this expression by using the elliptic theta function, 
and doing so, we reproduce the results obtained in Ref. 11). 

b) uj q 2 + a 2 - n 2 = (ttT) 2 + a 2 - /u 2 < 

In this case, the integrand in ()15|) is exponentially suppressed in the limit 
r — > —ioo for n > [N] and in the limit r — > +ioo for n < [N], where N(a) = 
(P^{i 2 -a 2 /Tr - l)/2 and [N] is the maximum integer smaller than N. Thus, we 
must add the residues of the poles for n < [N]. This yields 



I(a)+I*(a) 



1 



All 3 / 2 (3 
[N] 



:x - P OO 

V / dTf(a;n;T) 



n>[N] 



1/A 2 

(18) 



+ ho(a;n) +hj(a;n) + / d-r c/(ct; n; t)1 
<?(<r; n;r) = ^ coth(QHr) sm(2uj nf iT)e T ^ 2 -^ + ^ , (19) 

— J7r/4 /" 7r /2 

/j ( CT; „) = Qtfyl / e~ i6/2 cot(QHe id A 2 

2 J-tt/2 



7r/ < c<>m j i ii , 1 

-7T/2 

x exp[-i{(u; n - i//) 2 + o- 2 }e^ 2 ] + (c.c.) , (20) 

3/2 



fyto n ) = JL (^j e - 2 ^M/(Q^) 



n i 



xsin{^K 2 -^ 2 + f x 2 ) + ^ j . (21) 
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By choosing the contour of the proper time integration carefully, we find explicit 
expressions for the effective potential with Eqs. (16)-(21), without any ambiguities, 
even for fi > T > 0. This result is consistent with the previous results in the 
simultaneous limits H — ► and 1/A 2 — > 0. For example, the effective potential (|14|) 
exactly reproduces the result given in Ref. 5) at finite T and 



§3. Phase structure 



In this section, we plot the phase structure of the NJL model at finite T, \i and 
H. For this purpose, we numerically calculated the effective potential 1)14(1 . The 
parameters of our model are the coupling constant G and the proper time cutoff 
A. We fix these parameters to the physical mass scale that reproduces the observed 
values. Here, it is assumed that the T, /j, and H dependences of G and A are not 
so strong. Below, we use the values G = 38.7GeV~ 2 and A = 0.864GeV obtained in 
the appendix. 




H = 0.275GeV 



o [GeV] 

(a) At T = 0.05 GeV, H = 0.18 GeV 2 




o [GeV] 

(b) At T = 0.05 GeV, H = 0.3 GeV 2 



1 1 


1 1 1 J 




v - 0.322/ / 








n = 0.285 / ' 


1 


p = 0.270GeV 



a [GeV] 



(c) At T = 0.01 GeV, H = 0.18 GeV 2 . 
Fig. 2. Behavior of the effective potential with T, \i and H fixed. 



*' To compare the effective potential ffll4B with that obtained in Ref. 5), we must rewrite the 
equation using the renormalized coupling constant G r , which is defined by the renormalization 



condition 



<ry_ 

da 2 



1 

G~ r 
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jlx [GeV] 



Fig. 3. Dynamical fermion mass (a) as a function of the chemical potential fi at T = 0.01 GeV 
with H fixed at the typical values 0.0, 0.10, 0.18, 0.21 and 0.60 GeV 2 . 



Typical behavior of the effective potential is plotted in Fig. |2j As seen in the 
figures, we observe the followings: 

(a) There is a second order phase transition as [i is increased, with T kept large and 
H kept small. 

(b) There is a first order phase transition as \i is increased, with H kept large and 
fixed. 

(c) There are two steps of the transition as fi is increased, with T and H kept small. *) 
First, a first order transition takes place from the third non- vanishing extremum to 
the first extremum. Next, a second order phase transition occurs where the first 
non-vanishing extremum disappears. 

The ground state is determined by finding the minimum of the effective potential. 
The necessary condition for this minimum is given by the gap equation 



dV 



= 0, (22) 

a=(a) 



where (a) is the expectation value of a, which corresponds to the dynamically gen- 
erated fermion mass. If the gap equation (122|) has more than one solution, we 
numerically calculate the effective potential for each of these solutions and thereby 



*' Such a two-step transition is also found in the non-Abelian gauge theory in the Randall- 
Sundrum background. 22 ' 
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find the value (a) in the ground state. In Fig. El we plot the behavior of the dynam- 
ical fermion mass (a). We find that (a) disappears above a certain critical chemical 
potential. We clearly observe that (a) increases continuously as H increases for a 
small chemical potential. The magnetic field enhances the chiral symmetry breaking. 
This phenomenon is known as magnetic catalysis. 7 ^' 8 ^ For large //, the behavior of 
the dynamical fermion mass is different from that in the case of magnetic cataly- 
sis. Near the critical chemical potential, a variety of phenomena are observed. For 
H = 0.10 GeV 2 , there is an oscillation of the critical line. This is called the de 
Haas- van Alphen effect. 12 )' 23 ) It is caused when the Landau levels pass the quark 
Fermi surface. If H is large enough, the critical chemical potential increases as H 
increases, and a mass gap exists at the critical value of u. However, there is a com- 
plex behavior of the critical chemical potential for sufficiently small H. We observe 
the two steps of the transition at H = 0.18 GeV 2 . 

To understand the situations more precisely, we numerically calculated the crit- 
ical values of T, u and H at which the dynamically generated fermion mass disap- 
pears. For the second order phase transition, the critical point is obtained by taking 
the following limit of the gap equation: 

lim { — }=0. 

For the first order phase transition, we directly observe the behavior of the effective 
potential and find the critical point. We plot the phase boundaries in Fig. ^ On 
the T-u plane, it is clearly observed that the broken phase spreads out as H is 
increased for a small chemical potential. The first order phase transition disappears 
near H = 0.18 GeV 2 . 

We illustrate the details of the phase boundaries at large chemical potential in 
the H-T and fi-H planes. It is seen that the first and second order phase transitions 
coexist in the planes. As is seen in Fig. EJh), the broken phase is separated into 
two parts for u > 0.28 GeV. The distortion of the critical line at u = 0.34 GeV is 
caused by the de Haas- van Alphen effect. Such an effect is observed as an oscillating 
mode for T < 0.03 GeV in the u-H plane in Fig. Etc). In the parameter range 
corresponding to the second order phase transition, chiral symmetry breaking is 
suppressed as H increases, while it is enhanced in the range corresponding to the 
first order phase transition. 

§4. Summary 

We investigated the phase structure of the NJL model at finite temperature and 
chemical potential in an external magnetic field. The Fock-Schwinger proper time 
method was applied to a thermal field theory. Using the two-point function obtained 
in Ref. 13), we calculated an effective potential that exactly includes the effects of 
temperature, chemical potential and constant magnetic field to leading order in the 
1/N C expansion. We determined the physical contour of the proper time integral. 
There is no ambiguity with regard to the contour of the proper time integration for 
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(a) The T-p plot with H fixed. 



0.05 0.1 0.15 0.2 
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(b) The H-T plot with fixed. 
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(c) The \x- H plot with T fixed. 

Fig. 4. Critical lines in the T-/j plane, H-T plane and fi-H plane. The dashed curve represents the 
first order phase transition, while the solid curve represents the second order phase transition. 



[i > T > in the expression of the effective potential. 

We used the current algebra relationship with the current quark mass to de- 
termine the mass scale of the model. The position of the tricritical point depends 
strongly on the mass scale and the regularization procedure. Here, we used the 
parameter values 6? = 38.7 GeV -2 and A = 0.864 GeV, where the tricritical point 
appears in the T-\x plane. 

Through the study of the shape of the effective potential, we observed a phase 
transition from the broken phase to the symmetric phase when the temperature 
T, chemical potential [i and magnetic field H vary. Increasing the temperature 
and chemical potential tends to prevent the chiral symmetry breaking at H = 0, 
while increasing the magnetic field enhances the symmetry breaking at T = fj, = 0. 
Therefore, complex behavior of the phase boundary is found with the combined 
effects of T, fj, and H, as shown in Fig. 3. We found a two-step transition from 
the broken phase to the symmetric phase for a medium size magnetic field, H ~ 
0.18 GeV 2 . We found a separation of the broken phase into two parts in the H- 
T plane. With this separation, two tricritical points appear in the critical line for 
[i ~ 0.27 GeV. More than two tricritical points were observed in the critical line for 
T ~ 0.01 GeV in the [i-H plane. In the phase diagram, we observed a contribution 
from the magnetic field, which is known as magnetic catalysis, and the de Haas-van 
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Alphen effect. 

In the present paper, we regard a simple four-fermion interaction model as a low 
energy effective theory that may possess some fundamental properties of the chiral 
symmetry breaking in QCD. Although the present work is restricted to an analysis 
of the phase structure of the four-fermion interaction model with massless quarks, 
we hope to apply our result to critical phenomena in the real world. In QCD, the 
current quark mass for up and down quarks is much smaller than the constituent 
quark mass a for T = [i = H = 0. Because the expectation value of a decreases 
near the critical point, it is no longer valid to ignore the current quark mass. If the 
contribution from the current quark mass is included in the effective potential, the 
sharp transitions illustrated in Fig. 4 disappear, as shown in Refs.7) and 20). In that 
case, the second order phase transition becomes a cross over. 

The strength of the magnetic field that can affect the phase structure is larger 
than that existing in the core of neuron stars and magnetars. However, an oscillating 
mode has been observed even for small magnetic fields, and this may have some 
interesting contribution to the behavior of such stars. Furthermore, we believe that 
the magnetic field may have been significantly stronger in the early stages of our 
universe with the primordial magnetic field and that this may have contributed to 
the critical phenomena in the early universe. To investigate the phenomena exhibited 
by dense stars, we cannot avoid considering the vortex configurations that generate 
non-local ground state in the super-fluidity phase. It is necessary to calculate the 
effective action at finite [i and H. We will continue our work further and extend our 
analysis to non-local ground states. 
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Appendix A 

Physical Scale of G and A 

Here we discuss the physical scale of the parameters G and A of our model. 
We regard the NJL model as a low energy effective theory of QCD. The parameters 
should be determined so as to reproduce the observed quantities at low energy scale. 

First, we calculate the pion decay constant using the proper time method. The 
two-point function for a pion composite field is given by 



xtv[(-i 1 5 T a )S(x,0;k )(-i 7 5 r b )S(x',0-,k +p )}, (A-l) 
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where S(x, y; ko) is the quark Green function defined by 

(7% + il j dj - a)S(x, y; k ) = 5 3 (x - y). 



(A-2) 



Here, a is determined by finding the stationary point of the effective potential which 
is determined as a function of G and A. Using the proper time method, the Green 
function S(x,y;ko) can be written 



S(x,y;k ) 



e 7rj/4 

(4vr)3/2 



(iyidj + 7% + a) 



dr 



exp 



i (x-y ) 2 
4 r 



+ i(a 2 - k 2 )r 



,. (A ' 3) 

where r is the proper time. Substituting Eq. (|A-3|) into Eq. (jA- 1|) . and performing 
the integration over k, x and x' , we obtain 



No Nc_ 

G 27T 2 



+0 



(r + s) ; 



T + S 



2i-p* 



TS 



T + S 



a 



(A-4) 



The renormalization constant for the pion wave function is defined by r^\p 2 ) = 
Z~ l p 2 + 0(p 4 ). After applying the Wick rotation to Eq. (|A-4(I the renormalization 
constant is derived as 



2tt 



TS 



(r + s) 3 \r + s 



(A-5) 



To obtain the relationship between f n and Z^, we evaluate the order parameter 
of the chiral symmetry breaking, [iQg,7r 6 ], where Q1 is the conserved axial charge 
and 7r a is the asymptotic pion field. Then, using the conserved axial current j^(x), 
we obtain 







[iQiy] 0) = (0|[/d a xj o a (x),7r b ]|0) 



i3 =t -a 1 



(fx id^ ( 



The pion decay constant f n is defined by 

(0\f,(x)\ir(p)) = -if 7T 5 ab p fl e-^. 

This implies that 

Substituting Eq. ()A-8|) into Eq. (|A-6|) . we obtain 



o) = f^ ab . 



On the other hand, j^(x) and ir a (x) are described by the quark field ip(x) as 



G 



(A-6) 

(A-7) 
(A-8) 

(A-9) 
(A-10) 
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Fig. 5. Relationship between G and A 



Inserting Eq. (jA-10}) into the second line of Eq. HA-6|) , we find 



[iQlA 



0) s 



■ab 



-Z-V 2 a6 ab . 



(ATI) 



Comparing ()A-9|) with (|A-11|) . we obtain the relationship /„. = —Z n 1 ^ 2 c 
Eq. (|A-5|) and this relationship, it is found that the pion decay constant can be 
expressed as a function of o and A, 

U 2 = Q 2 dr± {e~^(l - a\) - A 2 Ei(Vr)} , (A-12) 

where Ei(— x) is the exponential-integral function. 

We determine the relationship between G and A such that the measured value, 
/tt = 93 MeV, is reproduced, as shown in Fig. [5J The phase structure in the T- 
[i plane differs dramatically for values of G and A on the solid curve (case i) and 
the dashed curve (case ii) in Fig. [21 Evaluating the effective potential, we plot 
the typical behavior of the critical line in T-fi plane in both cases. As shown in 
Fig. 113 only the second order phase transition takes place in case i, while the first 
and second order phase transitions coexist in case ii, and a tricritical point appears 
on the critical line. These results agree with the well-known result obtained using 
the other regularization procedure. 2 )' 5 )' 21 ) Therefore, we choose the mass scale in a 
range corresponding to case ii. In chiral limit, we determine only the curve in the 
G-A plane. 

To determine the values of G and A in this range, we evaluate the pion mass. 
As is well known, the pion is a massless Goldstone mode if the explicit chiral sym- 
metry breaking is not accounted for. In QCD, a small quark mass breaks the chiral 
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symmetry. The pion decay constant is proportional to a. Because the value of a 
is much larger than the current quark mass, we ignore the a contribution from the 
current quark mass when calculating the pion decay constant. However, to evaluate 
the realistic pion mass, m n = 138 MeV, we consider the influence of the current 
quark mass. 

Accounting for the current quark mass, the axial current is not conserved, and 
we have 

gPfflx) = mi>ij 5 T a 'ip(x), (A- 13) 

where m is the current quark mass, 2m = m u + m^. Substituting Eq. I)A-7|) into 
Eq. (|A-13[) . we obtain 

m (o |^7 5 r>(0)| 7r b (p = 0)) = -Uml5 ab . (A44) 
Similarly to Eqs. (|A-6|) and (|A-9|I . we obtain 

-U (O |^75T>(0)| Tr b (p = 0)) = (0 |^| 0> (A-15) 
in the chiral limit. Therefore, the pion mass is obtained as 

(/^) 2 = -^m(a), (A-16) 
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to first order in the current quark mass. 20 ) This is known as the Gell-Mann— Oakes— Renner 
relation. 

From the measured values f w = 93 MeV and m n = 138 MeV and Eqs. (|A- 12|) 
and (|A-16|) . we obtain the matching condition for the coupling constant G and the 
proper time cutoff A with the physical scale. These parameters are extracted as 
functions of rh. For rh = 5.5 MeV, the parameters G and A are found to be 

G = 25.4 GeV" 2 , A = 0.991 GeV, (A- 17) 

which are located in the range of the solid curve in Fig. \5\ A larger current quark 
mass is realized on the dashed curve in Fig. |5J If we set rh = 7.0 MeV in Eq. (lA~16|) . 
the parameters G and A are fixed as 

G = 38.7 GeV~ 2 , A = 0.864 GeV, (A-18) 

which are in the range of the dashed curve in Fig. |5J 

It should be noted that we consider the contribution from the current quark 
mass only for the purpose of determining the values of G and A. 
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